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In a recent letter [T], Amini et al. claim having found a 
mobility edge in graphene, a truly two-dimensional (2D) 
system. Their mobility edge ought to be induced by on- 
site uncorrelated (Anderson-type) disorder and — unlike to 
3D systems — shall separate localized states in the band 
center, from the remaining extended states. In order to 
distinguish localized from extended states, the authors an- 
alyzed the 'typical' density of states (DOS) [2], 



ftyp(^) = exp 
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(1) 



by means of the so-called 'Regularized Kernel Polynomial 
Method' (RKPM) [g. 

In this comment, we argue that the main conclusion 
of the above letter is wrong. Employing the Chebyshcv- 
based Kernel Polynomial Method (KPM) [i], we demon- 
strate that the Anderson model on the graphene lattice, 
with uniformly distributed on-site energies in the inter- 
val [— W/2, W/2], shows localization throughout the whole 
spectrum, even for small disorder strengths. That is, no 
mobility edge exists at all. 

In particular, we want to make the following points: (i) 
The RKPM is not well suited for the evaluation of the typ- 
ical DOS, addressing the subtle 2D Anderson localization 
phenomenon, (ii) A careful finite-size scaling is necessary 
to distinguish localized from extended states, either by 
means of the typical DOS or by considering the whole dis- 
tribution of the local density of states (LDOS), pi ~ p^. 

Let us substantiate these assertions in more detail: 

(i) The RKPM is a polynomial expansion based 
on Gaussian- broadened Legendre polynomials (Pn{x))a- 
This broadening is supposed to damp the Gibbs oscil- 
lations known from finite order polynomial expansions. 
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Fig. 1: (Color on-line) Upper panel: Comparison of pure ((7 = 
- dashed line) and broadened {a — 4/M - solid red line) 
expansion of a (5-function by Legendre polynomials of order 
M — 2000. The two insets show a magnification of the broad- 
ened expansion over different ranges. Lower panel: Minimum 
value of {5{x — a))a- for various parameters of a, M and u. 



However, the kernel of this method is not strictly positive 
(as opposed to the Chebyshev- based KPM with Jackson 
kernel |3]), i.e., the approximation of a positive function 
may become negative in some regions. This is illustrated 
in figUl which gives the Legendre polynomial expansion of 
a delta- function S{x ~ a) (upper main panel). The insets 
display magnifications of the broadened expansion, prov- 
ing the existence of negative values. To corroborate this 
fact we show by the lower panel of fig.[I]that the expansion 
oi d{x — a) contains negative values for various peak posi- 
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Fig. 2: (Color on-line) Left panel: Comparison of our KPM 
results for ptyp (pmo) with RKPM data taken from [T], using 
identical system size A'^ = 1000^, disorder strength W/t — 2.2 
and kernel width a — 8x 10''. Right panel: Finite-size scaling of 
normalized typical DOS ptyp/pme at constant resolution {Nk = 
500 states within the kernel). This resolution corresponds to 
the one of the RKPM data at the maximum of the DOS. The 
total number of realizations are KrXKs = 61440, 25600, 11712 
for N = 252^ 500^ 1000^ respectively 



tions a, expansion orders M and Gaussian broadening a. 
Although there is an optimal value crM, negative values 
persist for all aM . In addition, the underlying recurrence 
for the expansion coefficients becomes numerically insta- 
ble for a > 8/M and the expansion characteristics strongly 
depends on the function expanded. 

The strict positivity of the kernel is of fundamental im- 
portance for studying localization by means of the LDOS 
distribution approach. Artificial negative values of ps (E) 
will suppress the typical DOS ([T]) at energies where it 
would be normally finite. This issue causes the region 
of vanishing ptyp in the vicinity of the band center for 
W/t > 1.4 in fig. 3 of [l]. To pinpoint this drawback, in 
the left panel of fig. [21 we contrast the data of Amini et 
al. [I] with our results obtained by KPM with the posi- 
tive Jackson kernel using matching parameters. Our data 
do not show any vanishing ptyp- Obviously , Amini et al. 
wrongly attribute physics to a numerical deficiency. At 
this point, we like to emphasize that even for a correct 
positive kernel, ptyp may erroneously vanish if the kernel 
width becomes comparable to the level spacing. Since in 
graphene the mean level spacing strongly varies through- 
out the spectrum, an adaption of the kernel width to the 
mean DOS is necessary [5]. This guarantees that for any 
fixed energy at least some energy levels are within the 
width of the kernel for each realization of disorder. 

(ii) The use of only a single finite system size to dis- 
tinguish localized from extended states by means of ptyp 
is insufficient. In particular, a finite value of ptyp does 
not guarantee extended states. Instead, the scaling be- 
havior of the typical DOS should be considered (see fig. [21 
right-hand panel). If ptyp{E) / p^dE) remains constant 
for different system sizes, the corresponding states are ex- 
tended, whereas ptyp > indicates Anderson local- 
ization. According to our results all states in graphene are 




Fig. 3: (Color on-line) Finite-size scaling of the LDOS distri- 
bution with corresponding log-normal fit (left-hand panel) and 
normalized typical DOS ptyp/pmo (inset right-hand panel) at 
energy E/t = 0.25 for weak disorder strength W/t — 0.3. The 
total number of realizations are 4 x 10* for A'' = 1500^ and 10^ 
otherwise. 



localized for any disorder strength, in agreement with the 
one-parameter scaling theory of localization for 2D lat- 
tices [5]. Note that the sensitivity of ptyp is limited: for 
localization lengths much larger than the system size, a 
finite-size scaling of the whole LDOS distribution permits 
a more reliable detection of the localization properties [5] , 
offering two main advantages. First, reducing the distri- 
bution characteristics to ptyp rnakes only sense if the shape 
is close to log-normal. In more complex systems, this ap- 
proximation might be poor. Second, prior to a pronounced 
shift of the distribution maximum, i.e. a variation of ptyp, 
localization is indicated by a broadening of the LDOS dis- 
tribution with increasing system size (sec fig. [3]). As a 
consequence, our data for the LDOS distribution clearly 
indicate localization, also for weak disorder. 

To conclude, we have shown that under uncorrelated on- 
site disorder, graphene behaves like ordinary 2D systems, 
which show Anderson localization of all single-electron 
states. Hence, in contradiction with the findings by Amini 
et al. |T| , Dirac fcrmions will not remain delocalized below 
a critical disorder strength. The correct analysis of the 
typical DOS, supplemented by a careful finite-size scaling, 
is unsuggcstive of mobility edges around the Dirac point 
of graphene. 
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